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Abstract 

With 58 million produced J/ip events collected by the BES-II detector at the BEPC, the decays 
J/i/j — > AA and XrX are analysed. The branching ratios are measured to be Br(J/ijj — > AA) = 
(2.05 ± 0.03 ± 0.11) x 10~ 3 and Br(J/tp -► = (1.40 ± 0.03 ± 0.07) x 10~ 3 . The angular 

distribution is of the form = 7V (1 + acos 2 0), with a = 0.65 ± 0.12 ± 0.08 for J/i/> -► AA and 

a = -0.22 ± 0.17 ± 0.09 for J/ip -► respectively. 

PACS numbers: 13.20.Gv, 14.20.Jn, 23.20.En 
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As is well-known from the helicity formalism, the angular distribution of B in the decay 
of a neutral vector resonance V into a baryon-antibaryon pair BB is given by Q 

dN 



dcosO 



1 + acos 9, 



where 9 is the emission polar angle of B in the V rest frame. The first order calcula- 
tions 

as 

of the perturbative QCD predict the theoretical value of a at J Jib energy. Table 
1 summarizes the theoretical predictions of a for the decays J lib — > AA, X°£ . Several 
experiments have measured a for J Jib — > pp Q Q Q, AA,E°E° Q and E~E + fl. 

In previous papers Q 3, the analyses of J Jib — > AA and E°E° using 7.8 x 10 6 J/ip events 
collected with the BES-I detector have been reported. The BES-I value of AA angular dis- 
tribution coefficient a is in good agreement with the DM2 value [7] as well as the theoretical 
predictions in Ref.[3]. However the BES-I value of angular distribution coefficient a 

with minus sign is obviously deviated from those of DM2 , MARKII as well as the theoretical 
prediction. 

This letter presents the analyses of the decays J Jib — > AA and based on BES-II 

58 x 10 6 J Jib events, with the purpose of improving the accuracy of branching ratio and a 
value measurements and clarifying the sign of the a for J Jib — > 

Table 1 theoretical expectations for am J Jib BB. 



Channel 


Ref.[2] 


Ref.[3] 


AA 


0.32 


0.51 


S°E° 


0.31 


0.43 



The BES-II detector has been described in detail elsewhere jlo| . 

Since the decays studied include AA pair in the final state, selection criteria are 
used to select J Jib — > AA + X events, where X is a system of neutral particle(s) and/or 
undetected charged particle(s). The A is identified by its pir~ decay mode. Candidates 
for J Jib — > AA + X events are selected by requiring exactly four reconstructed charged 
tracks in the drift chamber with zero net charge. Tracks with | cos9 c h |< 0.8 and transverse 
momemtum p xy > 0.07 GeV are accepted, where 9 c h is the polar angle with respect to the 
beam direction. Because the A is produced at the second vertex, no limits are required for 
the primary vertex position. The particles are identified by requiring that their combination 
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weights of the time-of-flight (TOF) and the ionization energy loss (dE/dx) in the drift 
chamber be consistent with the corresponding particle hypothesis. 

The events with four charged particles satisfying hypothesis of a pair of pp and a 
pair of 7r + 7r~ are selected. Fig. 1 shows pn~ invariant mass distribution for remaining 
J/ip — > pir~piT + + X candidates (solid line). By fitting the M pw - distribution to a gaussian 
distribution plus a quadratic polynomial, a 2.6 MeV mass resolution is estimated for the 
peak. Fig. 1 also shows pn~ invariant mass for Monte Carlo simulation (dashed line) with the 
mass resolution of 2.6 MeV. To select A and A, the requirements of | M p7T - — 1.1156 |< 0.008 
GeV and | M p ^+ - 1.1156 |< 0.008 GeV are imposed. 

The energy distribution of the AA pair, E AA - ,for J/ip — > AA + X candidates is shown 
in Fig. 2, where E AA is defined as the energy sum of A pair, data is represented by solid line 
and Monte Carlo simulation by dashed line. 

The contamination from non-AA+A" events in Fig. 2 can be estimated by the contribution 
from the sidebands, namely the off- A resonance background shown in Fig. 1. The events 
in the sidebands are defined as those in the squre area of | M im - - 1.1156 |< V2 x 3a 
and | M p7r + — 1.1156 |< V2 x 3a minus the squre area of | M p7T - — 1.1156 |< 3a and 
| M pn +_ — 1.1156 |< 3a in two dimentional M pn - vs. M pn + plot, where a = 0.0026 GeV is 
the mass resolution of A (A). In Fig. 2 the sideband contribution is indicated by shaded 
area. 
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Fig. 1. pn~ invariant mass distribution in J/ip — > p7r~pir + + X candidates for data 
(solid line). The dashed line is Monte Carlo simulation for J/ip — > A A. 
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Fig. 2. Energy distribution of AA for J/ip —> AA + X candidates. The shaded area is 
the sideband contamination described in the text. The dashed lines are from Monte Carlo 
simulation of J/ip — > S°E° and J/ip —> AA. 




A clear peak centered at the J/ip mass in Fig. 2 from the decay J/ip — > AA is observed. 
The enhancement centered at 2.9 GeV is due to the decay J/ip — > £°£°, where S° — > A7. 
To obtain J/ip ^ AA from the selected J/ip —> AA + X candidates, contamination from 
J/ip decay is suppressed by imposing the total missing momentum cut p m j SS < 0.15 

GeV and the A (A) momentum requirement 0.98 < p A < 1.2 GeV. Finally by requiring 
3.02 < E Kli < 3.2 GeV a AA sample of 9462-503=8959 events is obtained, here 503 events 
are from the sidebands, which are obtained using the same criteria for AA selection to the 
aforementioned non-AA + X sideband events. 

The detection efficiency for the events J/ip —> AA — > p7r~p7r + depends on the A direction. 
A Monte Carlo simulation gives the detection efficiency e{6j) in different A polar angle 9i 
with the bin size Acosdi = 0.1. The efficiency corrected angular distribution of A for the 
decay J/ip — > AA is shown as histogram line in Fig. 3. Fitting this angular distribution to 
the theoretical form 

N (l + acos 2 6), 



dcosO 
yields 

N = (1990.9 ± 40. 1)/0.1, 
a = 0.65 ±0.12, 

where the error is statistical. 
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Using the Nq and a values, the number of J /if} — > AA events corrected for the efficiency 
6j is deduced from the equation N AA = A r (l + acos 2 9)dcos9, and the branching ratio is 
obtained by the equation 



5r(J/V> -> AA) 



N. 



AA 



Br (A — > p7t~)Nj/^ 



where A r j/^ ) = 58 x 10 6 is the number of J/ip events with the error of 4.72% 11 1. 

The remained contamination from J/if> — > is 0.8% using the branching ratio deter- 
mined by this work. Subtracting this contamination the number of J /if} — > AA signal events 
is 8887. The branching ratio and angular distribution coefficient are obtained to be 

Br(J/ip -> AA) = (2.05 ± 0.03 ± 0.11) x 10~ 3 , 



a 



0.65 ±0.12 ±0.08, 



where the first error is statistical and the second error is systematic. The systematic error 
includes the uncertainty of the detection efficiency due to imperfection of Monte Carlo 
simulation, the contamination from the indefinite branching ratio of background channel 
J /if) — > S°A ± c.c. and the uncertainty of the number of J/tf> events. 
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Fig. 3. The A angular distribution for J/if> —>■ AA events. The histogram is efficiency 
corrected data. The smooth dashed line is the fitting result. 



The enhancement centered at 2.9 GeV in Fig. 2 is dominantly due to the decay J /if} — > 
E°E°,where S° — > A7. 2.75 < E A -^ < 3.02 is required to select E°E° events. To remove the 
background from J/ip — > AA, AA7r°, AA7 and S°A± c.c, the events are kinematically 
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fitted to the J/ip — > AA77 topology by imposing energy and momentum constraints (4C). 
The combination with the smallest x 2 i n the 4C fit is chosen to identify the radiative photons 
if there are more than two photon candidates in an event, here a photon is defined as a 
cluster with deposite energy larger than 30 MeV in the barrel shower counter, outside a 
25° cone around p and outside a 12° cone around the other charged particles. Then a 6c 
fit ( two additional constraints Ma 7 = M^o, = M^o) is performed. It is required that 
^aa 77 (^ c ) < ^AA7^ C )' -^a 7 distribution is shown in Fig. 4, where the solid line is data 
and the dashed line is Monte Carlo simultaion J/ip — > The Ma 7 distribution is similar 

to Fig. 4. Finally, after applying X A a 7T (6c) < 40, 2194-82=2112 events are remained, where 



1^0 



82 is from the sidebands which are obtained using the same criteria for E°E" selection to 
the aforementioned non-AA + X sideband events. 



Table 2 The background for J/ip — > S°E 



Channel 


contamination (%) 


AA 


0.1 




0.3 


AAtt 


0.5 


E°A + c.c. 


< 0.08 


AA7 


< 0.08 



The detection efficiency for the signal events depends on the S° direction. The angular 
dependence of S°E° detection efficiency is obtained by selecting Monte Carlo events J/ip — > 
— > A7A7 — > j97r _ 7p7r + 7 through the same selection criteria. The efficiency corrected 
angular distribution is shown in Fig. 5 as histogram. Fitting this angular distribution to 
the theoretical form 

= N 1 + acos 2 e , 

dcosO 

yields 



AT = (1786.9 ±57.4)/0.1, 
-0.22 ±0.17, 



a 



where 9 is the polar angle between the S° and e + beam and the error is statistical. The 
number of J/ip — > events corrected for the detection efficiency is A^^o = N (l + 



8 



acos 2 8)dcos8 and the branching ratio is calculated with 

Br(J/ip 



Br 2 (A^pn-)Nj^ 



Table 2 summarizes the remaining contamination to J/ip — > S°S candidates from back- 
ground channels. The branching ratio of J ftp — > AA is taken from this work, while the 
others are from the PDG. Subtracting all the contaminations of first three channels, the 
number of J/tp — > signal events is 2093. The backgrounds of other two channels are 

considered in the systematic error. The results for the branching ratio of J/tj) — > S°S° and 
its angular distribution coefficeincy are 



Br(J/ip -> £°I7) = (1.40 ± 0.03 ± 0.07) x 10 



-3 



a = -0.22 ±0.17 ±0.09, 

respectively, where the first error is statistical and the second systematic. The systematic 
error includes the uncertainty of detection efficiency due to imperfection of Monte Carlo 
simulation, the contaminations from the indefinite branching ratios of background channels 
J/ip — > AA7 and J/ip — > S°A ± c.c. and the uncertainty of the number of J/ip events. 

Table 3 and Table 4 summarize the branching ratios and a values measured by this 
work for the decays J/ip — > AA and J/tp — > together with those previously reported 

by Mark II Q and DM2 Q. 
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Fig. 4. The S signal after 4c fit X^(4c) < 40. The solid line is data and the dashed line 



is Monte Carlo simultaion J/ij) — > S°S . 
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Fig. 5. The S° angular distribution for J/ip — > The histogram is efficiency corrected 

data. The smooth dashed line is the fitting result. 



The statistics of the signal events in BES-II result is much better than those of Mark II, 
DM2 and BES-I. In this experiment the a value for J /if} — > AA obtained by this experiment 
is consistent with the previous measurements, while the a value is negative for J/ip —> S°S°, 
which agrees with that of BES-I P| while conflicts with the values of MARK-II and DM2 
and theoretical expectation 3] Q 

In this analysis the 6c kinematic fit is used to select J lip —> £°£° and hence the contam- 
ination is small; while in DM2 p| and BES-I analyses [9j the S° is not reconstructed and 
cosO of S° is replaced with that of A. This makes a for J/ip — > S°E° determined by this 
work rather believable. 

It is worth noting that the central value of a reported by DM2 7] for J/if> —>■ is 
positive, however, the angular distribution has a convex shape as a whole, with the only 
exception at cos 9 = —0.65. 



Table 3 



Experimental measurements for decay J/ip — > AA 



Exp. 


a 


5r(xl0- 3 ) 


Evts 


MARKII 


0.72 ±0.36 


1.58 ±0.08 ±0.19 


365 


DM2 


0.62 ±0.22 


1.38 ±0.05 ±0.20 


1847 


BES-II 


0.65 ±0.14 


2.05 ±0.03 ±0.11 


8887 
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Table 4 



Experimental measurements for J/ip — > S°S channel 



Exp. 


a 


5r(xl0- 3 ) 


Evts 


MARK- 1 1 


0.70 ± 1.10 


1.58 ±0.16 ±0.25 


90 


DM2 


0.22 ±0.31 


1.06 ±0.04 ±0.23 


884 


BES-II 


-0.22 ±0.19 


1.40 ±0.03 ±0.07 


2093 
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